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METHOD FOR FILTERING SEISMIC DATA, PARTICULA RIX-BX, 
KRIGING _ 

GENERAL TECHNICAL FIELD - BACKGROUND ON KRIGING ANALYSIS 

5 The present invention relates to filtering seismic 

data, in particular by kriging analysis. 

Kriging analysis enables a random function to be 
resolved from its covariance function. 

In particular, it is conventionally used in 
10 geostatistics for filtering seismic data, particularly 
but in non-limiting manner, for characterizing 
reservoirs. 

Kriging analysis relies in particular on the 
assumption that a phenomenon measured locally by means of 

15 optionally-regular sampling can be analyzed as a linear 
sum of a plurality of independent phenomena,- the 
variogram of the overall phenomenon corresponding to the 
linear sum of the variograms of each of the independent 
phenomena making it up. 

20 .Conventionally, the variogram corresponding to the- 

measured experimental data is resolved as sum of modeled 
variograms, and from the experimental data and the models 
selected for the individual variograms -used to resolve 
the data, there are deduced the individual functions that 

25 make up the random function corresponding to the overall 
phenomenon. 

It is thus possible to extract from a seismic data 
map of the type shown in Figure 1 (e.g. raw experimental 
data) firstly the white noise that is present in the data 

30 (Figure 2a), secondly noise corresponding to linear lines 
(Figure 2b), and finally filtered data cleared of both of 
these kinds of noise (Figure 2c) . ' 

The kriging calculations for determining the values 
of the individual functions into which an overall random 

35 function is resolved are themselves conventionally known 
to the person skilled in the art. 
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In this respect, reference can be made, for example, 
to articles and publications mentioned in the 
bibliography given at the end of the present description. 

Very generally, the value of an individual function 
involved in making up the overall random function is 
determined as being a linear combination of experimental 
values for points in an immediate neighborhood of the 
point under consideration, these experimental values 
being given weighting coefficients. 

In other words, if it is considered that a function 
Z (x) is made up as the sum of individual functions Y u (x), 
this can be written: 

and the component Y u (x) is estimated by: 

a=\ 

where a is a dummy index designating the points under 
consideration around the point x for which it is desired 
to determine the estimated value Y u *(x), Z x being the 
value at the point x, N being the number of such points. 

It can be shown that the weighting coefficients \<x 
satisfy the equation: 
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where the index 0 designates the point for which an 
estimate is to be determined, the values to Cq N being 

the covariance values calculated from the model u 
corresponding to the component Y u (values of the 
covariance function for the distances between each data 
point and the point to be estimated) , the values C ±j being 
covariance values calculated as a function of the 
selected model for the variogram of the function to be 
estimated (values of the covariance function for the 
distances between the data points) . 
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These weighting coefficients A, x u are thus determined 
merely by inverting the covariance matrices. 

PROBLEMS POSED BY THE STATE OF THE ART - SUMMARY OF THE 
5 INVENTION 

One of the difficulties of presently-known kriging 
analysis techniques is that they require the use of 
models of the covariance functions. 

The advantage of using such models is that they make 
10 ' it possible to have matrices which are defined, positive, 
and invertible. 

Nevertheless, it will be understood that although 
such filtering techniques give good results, they are 
strongly dependent on the individual expertise of the 
15 person selecting the models for the various variograms . 

That can be a source of error, and prevents those 
techniques being used by people who are not specialists. 

Furthermore, selecting models also leads to 
significant losses of time in production. 
20 An object of the invention is to mitigate that 

drawback and to propose a filtering technique using 
kriging analysis that can be implemented in automatic or . 
almost automatic manner. 

The invention provides a method of filtering at 
25 least two series of seismic data representative of the 

same zone, by determining (e.g. by determining the cross 
variogram of the data series and solving the co-kriging 
equation) an estimate of the component that is common to 
the data series, and deducing a resolution of these data 
30 . series is deduced from the estimate. 

The invention also provides a method of processing 
seismic data in which a filter method of the above- 
specified type is implemented in order to compare two 
series of seismic data corresponding, for the same zone, 
35 to grids of at least one common attribute obtained for 
two distinct values of at least one given parameter. 
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The invention also provides a method of filtering at 
least one series of data representative of the values of 
at least one physical parameter over at least one zone, 
characterized by identifying a model of a component of 
5 three-dimensional variability of its variogram, 

substracting said model from the experimental, variogram, 
and solving the kriging equation corresponding to the 
r different variograms in order to deduce an estimate of 
the corresponding variability component on the data 
10 series. 

BRIEF DESCRIPTION OF THE FIGURES 

Other characteristics and advantages of the 
invention appear further from the following description 
15 which is purely illustrative and non-limiting and should 
be read with reference to the accompanying drawings, in 
which: 

• Figures 1 and 2a, 2b, and 2c, described above, 
illustrate an example of seismic data mapping and of the 

20 corresponding resolution by kriging analysis; 

• Figures 3a and 3b show two maps of the same zone, 
obtained from acquisitions undertaken at two different 
times; 

Figure 4 is .a map of the component that is common 
25 to, the maps of Figures 3a and 3b; 

• Figures 5a and 5b and Figures 6a and 6b are maps 
of components other than the component that is common to 
the maps of Figures 3a and 3b; and 

• Figures 7a to 7c are graphs showing the 

30 distribution of errors respectively when using standard 
filtering, filtering by conventional factorial kriging, 
and filtering as proposed by the invention (factorial co- 
kriging) . 
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DESCRIPTION OF ONE OR MORE IMPLEMENTATIONS OF THE 
INVENTION 

Automatic filtering 

It is assumed that two maps are available that have 
5 been obtained for the zone with seismic data acquired, 
for example, at different instants or for seismic 
attributes that are different. 

By way of example, these two maps are of the type 
shown, in Figures 3a and 3b. 
10 The two functions corresponding to these two data 

series are written Zl and Z2 below. 

It is proposed to resolve each of these two . 
functions into the sum of their common component plus 
orthogonal residues . . 
15 For this purpose, there is determined, from two data 

series for which a cross variogram is available having 
the following values: 

1 v 

Yi 2 (h)=— X (Zl(x)-Zl(x+h))(Z2(x)-Z2(x+h)) 

where x and x+h designate the pairs of points taken into 
20 consideration in the direction and for the distance h for 

which the value of the variogram is determined, and where 

N is the number of pairs of points for said direction and 

said distance. 

Knowing this cross variogram, an estimate is then 
25 determined of the function corresponding thereto, which 

satisfies: 

where a and P are two dummy indices designating the 
points taken into consideration around the point x for 
30 which it is desired to determine an estimate of said 

function, Z l a and being the value at said point x./ N 

being the number of said points, and where X ] a and X 2 fi are 

weighting coefficients . 

These weighting coefficients X\ and X 2 fi are 

35 determined by inverting the co-kriging equation: 
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where the coefficients C12aP and C21ap are the cross- 
variance values of the functions Zl and Z2 at the points 
5 corresponding to indices a and P, and where the 

coefficients CllaP and C22aP are the covariance values 
respectively of the function Zl and of the function Z2 at 
said points. The index X corresponds to the point 
referred to above as x. 

10 It should be observed that the matrix which appears 

in this equation has the advantage of being invertible 
under certain calculation conditions. 

. In this way, using experimental covariances, the two 
variables corresponding to the two initial data series 

15 are resolved automatically into a common component and 

into two residual orthogonal components. The regularity 
of the data means that the experimental covariance is 
known for all of the distances used, so no interpolation 
is needed, so the matrix is defined positive. 

20 The function then obtained is an estimate of the - 

component that is common to both data series. 

Figure 4 shows an estimate of the common component 
obtained from the data corresponding to the maps of 
Figures 3a and 3b. 

25 It will be understood that it is particularly 

advantageous in numerous applications, and in particular 
in 4D seismic surveying, to have this common component 
available. 

It serves in particular to eliminate non-repeatable 
30 acquisition artifacts: 
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• on the basis of seismic attribute grids, and in 
particular, for example, on the basis of root mean square 
(rms) amplitudes in an interval; 

• on the basis of seismic time grids, and for 
5 example on the basis of seismic event times; and 

• on the basis of seismic velocity volumes, and for 
example on the basis of automatically-determined points 
of origin for velocity vectors. 

It may also be used in simple seismic surveying to 
10 eliminate non-repeatable acquisition artifacts: 

• in particular on the basis of seismic attribute 
grids calculated on consecutive incidence classes, or 
indeed 

• on the basis of seismic attribute grids calculated 
15 on volumes obtained by partial summing or converted 

waves . 

Furthermore, once an estimate has been determined 
for the common component, it is possible to determine the 
residual components corresponding to the difference 
20 . between the initial data and said estimated common 
component . 

These residual components can themselves be resolved 
by kriging analysis. 

This is shown in Figures 5a, 5b and 6a, 6b which are 
25 maps of white noise and of linear line noise estimated in 
this way for each of the two series of measurements shown 
in Figures 3a and 3b. 

Examples of error measurements as obtained by 
standard filtering, by kriging analysis filtering, and by 
30 co-kriging analysis filtering (or multivariable kriging) , 
are illustrated by the graphs of Figures 7a to 7c. 

On reading these figures, it will be understood that 
filtering by co-kriging analysis makes it possible to 
obtain dispersions that are much smaller than those 
35 obtained with conventional filtering or filtering by 
kriging analysis, and gives better results. 
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The above description relates to an implementation 
using two series of data. 

As will be readily understood, the proposed method 
can also be implemented in the same manner with a larger 
5 number of data series (campaigns) . 

Semi-automatic filtering 

This second implementation also enables simplified 
resolving when only one seismic data series is available 
10 (function SI) . 

It assumes that a model of a component 6m of the 
experimental variogram 0 is previously available. 

This model which is known beforehand is, for 
example, the model of an independent component of the 
15 subsoil geology: white noise, stripes, etc. 

Knowing this model of the component 0m, there is 
deduced therefrom the residual variogram corresponding to 
the difference between the experimental variogram and 
this component 0m. 
20 Kriging analysis is then performed in order to 

determine firstly the model component Sm and secondly, on 
the basis of the residual variogram, the orthogonal 
residue Rl such that : 

SI = Sm + Rl 

25 This automatic resolution . enables acquisition 

anomalies to be filtered out when they present three- 
dimensional coherence that is easily identified and 
modeled, such as the stripes parallel to the cables that 
are observed in the amplitudes and times when performing 

30 offshore seismic surveys. 
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